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SUMMARY 


The  dynamic  mechanical  properties  of  multi-constituent  particulate  composites,  consisting  of 
individual  Ni  and  A1  particles  dispersed  in  an  epoxy  matrix  are  investigated  in  this  study. 
Properties  of  such  composites  depend  on  the  mechanical  and  physical  properties  of  the 
individual  components;  their  loading  density;  the  shape  and  size  of  the  particles;  the  interfacial 
adhesion;  residual  stresses;  and  matrix  porosity.  These  multi-phase  particulate  composites 
systems,  particularly  those  with  high  fdl  densities,  have  not  typically  been  studied  rigorously,  to 
date.  Investigation  of  the  effects  of  higher-order  microstructural  features,  such  as  particulate  size, 
dispersion,  etc.,  on  the  static  and  dynamic  mechanical  response  of  these  multi-phase  (n>2) 
polymer-metal-composites  was  performed  using  a  factorial  design  of  experiments.  The  high 
strain  rate  compressive  properties  of  these  materials  were  characterized,  using  a  split  Hopkinson 
pressure  bar,  and  the  elastic  properties  of  these  complex  composites  were  determined  using 
dynamic  mechanical  analysis.  The  properties  were  correlated  with  microstructural 
characteristics  using  factorial  design  concepts  to  establish  the  effects  on  strength  at  low  and  high 
strain  rates. 

This  report  represents  the  completion  of  an  AFOSR  funded  effort  on  epoxy-based  multi-phase 
particulate  composites.  Several  published  conference  papers  and  presentations  follow  to 
summarize  the  results  of  this  effort. 
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ABSTRACT 

The  dynamic  mechanical  properties  of  multi-constituent  particulate  composites,  consisting  of 
individual  Ni  and  A1  particles  dispersed  in  an  epoxy  matrix  are  investigated  in  this 
study.  Properties  of  such  composites  depend  on  the  mechanical  and  physical  properties  of  the 
individual  components;  their  loading  density;  the  shape  and  size  of  the  particles;  the  interfacial 
adhesion;  residual  stresses;  and  matrix  porosity.  These  multi-phase  particulate  composites 
systems,  particularly  those  with  high  fill  densities,  have  not  typically  been  studied  rigorously,  to 
date.  Investigation  of  the  effects  of  higher-order  microstructural  features,  such  as  particulate  size, 
dispersion,  etc.,  on  the  static  and  dynamic  mechanical  response  of  these  multi-phase  (n>2) 
polymer-metal-composites  was  performed  using  a  factorial  design  of  experiments.  The  high 
strain  rate  compressive  properties  of  these  materials  were  characterized,  using  a  split  Hopkinson 
pressure  bar,  and  the  elastic  properties  of  these  complex  composites  were  determined  using 
dynamic  mechanical  analysis.  The  properties  were  correlated  with  microstructural  characteristics 
using  factorial  design  concepts  to  establish  the  effects  on  strength  at  low  and  high  strain  rates. 

INTRODUCTION 

Particulate  composite  materials  composed  of  one  or  more  varieties  of  particles  in  a  polymer 
binder  are  widely  used  in  military  and  civilian  applications.  They  can  be  tailored  for  desired 
mechanical  properties  with  appropriate  choices  of  materials,  particle  sizes  and  loading  densities. 
Several  studies  on  similar  epoxy-based  composites  have  been  reported  and  have  shown  that 
particle  size  [1,  2],  shape  [3],  and  concentration  [4]and  properties  of  the  constituents  can  affect 
the  properties  of  particulate  composites.  In  composites  of  AI2O3  particles  in  epoxy  (Epon 
828/Z),  increasing  the  particle  concentration  and  decreasing  the  particle  size  were  found  to 
increase  the  stress  at  4%  strain  [5].  A  study  of  aluminum  filled  epoxy  (DGEBA/MTHPA)  found 
adding  a  small  amount  of  filler  (~  5  vol.%)  increased  the  compressive  yield  stress,  but  additional 
amounts  of  filler  decreased  the  compressive  yield  stress  [6].  However,  tests  on  epoxy  (DOW 
DER  331/bisphenol-A)  found  that  increasing  the  volume  percent  of  glass  bead  filler  increased 
the  yield  stress  and  fracture  toughness  of  the  material  [7,  8].  In  a  study  on  a  similar  material, 
decreasing  the  aluminum  particle  size  from  micro  to  nano  resulted  in  increased  epoxy  crosslink 
density  and  subsequently  increased  static  and  dynamic  strength  [1].  In  this  study,  a  factorial 
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design  of  experiments  is  used  to  examine  the  effect  of  aluminum  particle  size  and  aluminum  and 
nickel  volume  percents  in  epoxy-based  Al-Ni  particulate  composites. 

Design  of  experiments  (DOE)  has  been  extensively  used  to  optimize  processes  in  a  wide  variety 
of  fields  [9].  However,  the  traditional  approach  in  mechanical  properties  testing  has  been  to 
change  one  factor  at  a  time,  which  can  be  expensive,  time  consuming,  and  does  not  reveal  the 
interactions  between  two  factors.  A  factorial  design  of  experiments,  in  which  all  possible 
combinations  of  the  levels  of  the  factors  are  investigate  in  each  replicate,  provides  an  efficient 
way  to  fully  understand  the  effects  of  the  individual  factors  as  well  as  their  interactions  with  the 
other  factors  [9].  The  simplest  DOE  is  a  2^  factorial  design,  where  each  factor  (1,  2,  ...,  k)  is 
given  a  low  (-1)  and  high  (-1-1)  level  [10],  where  the  levels  -1  and  -i-l  are  the  coded  factors.  With 
this  design,  a  linear  model  for  all  factors  and  their  interactions  can  be  developed.  More  complex 
designs  or  additional  levels,  such  as  centerpoints,  can  be  used  to  determine  the  nonlinear  or 
higher  order  effects.  A  2^  factorial  design  could  require  an  excessive  number  of  runs,  if  k 
becomes  very  large.  For  example,  a  two-level  factorial  design  with  6  factors  requires  64  runs. 
Depending  on  the  system,  a  fractional  factorial  may  be  appropriate,  where  only  selected  cases  are 
run  but,  by  careful  selection  of  those  cases,  main  effects  and  even  two  level  interactions  can  be 
investigated  [9,11]. 

The  output  of  a  2*^  design  of  experiments  is  a  model  for  the  system  with  the  low  and  high  levels 
of  all  of  the  factors.  For  a  k  =  2  experiment,  the  model  can  be  written  as 


y  =  /?o  +  /?i^i  +  /?2^2  +  1^12X1X2  +  e  (1) 

where  y  is  the  response,  xj  and  X2  are  the  coded  factors,  and  s  is  the  error  [A].  In  order  to 
determine  which  of  the  factors  and  interactions  are  significant,  an  analysis  of  variance  (ANOVA) 
is  used,  which  partitions  the  variability  into  its  component  parts,  i.e.  factors,  interactions,  or  error 
[9].  The  analysis  computes  an  F-value  and  a  probability  of  achieving  that  F-value  (p-value), 
which  is  how  likely  is  it  that  this  variability  due  to  a  particular  factor  is  due  only  to  noise.  A  high 
F-value  and  low  p-value  indicates  that  a  factor  is  significant  to  the  model.  From  the  models 
developed  with  DOE,  process  can  be  optimized  or  factors  that  have  no  effect  can  be  chosen 
based  on  other  considerations,  such  as  cost. 

EXPERIMENTAL  SET-UP 

Composites  of  aluminum  and  nickel  powders  in  an  epoxy  binder  were  prepared  for  this  study. 
The  aluminum  particle  size  was  varied  between  5  (Valimet,  H5  aluminum)  and  50  (Valimet,  H50 
aluminum)  pm.  The  H5  aluminum,  shown  in  Figure  1(a),  was  found  to  have  an  average  particle 
size  of  5.43  pm  with  spherical  smooth  particles.  The  H50  aluminum  particles,  Figure  1  (b),  were 
also  smooth  and  nominally  spherical  with  an  average  particle  size  of  51.91  pm.  The  nickel 
particles  (Micron  Metals)  had  rougher  surface  texture  and  more  irregular  shape  shown  in  Figure 
1  (c),  with  an  average  particle  size  of  47.45  pm.  However,  there  was  also  a  small  fraction  of 
particles  in  the  nickel  powder  that  had  an  average  particle  size  of  97.44  pm. 
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Figure  1.  SEM  micrographs  of  (a)  H5  aluminum,  (b)  H50  aluminum  and  (c)  nickel 
powders 


Table  1.  Factorial  design  of  experiments  for  epoxy-based  particulate  composites  containing 
aluminum  and  nickel 


Sample 

A1  Particle  Size  (pm) 

A1  Volume  Percent 

Ni  Volume  Percent 

MNML-8 

5 

20 

0 

MNML-4 

5 

20 

10 

MNML-6 

5 

40 

0 

MNML-2 

5 

40 

10 

MNML-7 

50 

20 

0 

MNML-3 

50 

20 

10 

MNML-5 

50 

40 

0 

MNML-1 

50 

40 

10 

The  composite  materials  used  in  this  study  were  prepared  using  a  factorial  design  of 
experiments,  with  the  variables  being  the  aluminum  particle  size  (5  or  50  pm),  the  volume 
percent  of  aluminum  (20  or  40  vol.%),  and  the  addition  of  nickel  (10  vol.%).  The  full  factorial 
design  of  experiements  is  presented  in  Table  1. 

Compression  experiments  at  quasi-static  strain  rates  were  conducted  with  an  Instron  5500  testing 
system  with  a  5  kN  test  frame.  Care  was  taken  to  center  the  samples  on  the  platens  prior  to 
testing.  Instron  Bluehill  software  was  used  to  conduct  constant  extension  tests  at  approximately 
4x10'"^  /s.  A  thin  film  molybdenum  disilicide  (MOSii)  was  used  to  lubricate  the  surfaces  of  the 
platen  in  contact  with  the  test  specimen.  The  strain  in  the  sample  was  measured  using  an  Instron 
Advanced  Video  Extensometer  (AVE),  in  which  the  movement  of  small  dots  placed  on  the 
sample  are  tracked. 

For  MNML-2  and  MNML-3  samples,  the  5  kN  load  frame  on  the  Instron  5500  was  sufficient  to 
load  the  materials  approximately  to  yield.  However,  at  this  point  the  strength  of  the  material 
overwhelmed  the  load  cell.  These  materials  were  tested  past  yield  at  the  same  strain  rates  using 
an  MTS  810  testing  system  with  a  100  KN  test  frame.  The  details  of  this  testing  are  described  in 
a  previous  paper  [12].  That  the  measured  stress-strain  curves  in  the  quasi-static  regime  are 
serrated  in  nature  is  believed  to  be  an  artifact  of  the  testing  set-up  rather  than  due  to  any  intrinsic 
material  property.  The  tests  were  conducted  under  constant  true  strain  conditions  by 
continuously  modulating  the  ram  speed  rate  according  to  the  measured  force  and  ram 
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displacement,  in  a  feedbaek  loop.  Faetors  such  as  machine  eompliance  and  stick-slip  loading  on 
the  compression  speeimen  faces,  along  with  any  visco-elastic  behavior  of  the  epoxy  matrix 
material  can  dramatically  affect  the  eontrol  feedbaek  to  the  ram,  resulting  in  serrated  loading 
curves. 

Compression  experiments  at  intermediate  strain  rates  (approximately  1x10^  and  5x10^  s'^)  were 
conducted  using  a  split  Hopkinson  pressure  bar  (SHPB)  [13]  system.  The  experiments  were 
conducted  using  the  SHPB  system  loeated  at  AFRL/RWME,  Eglin  AFB,  FL,  whieh  is  comprised 
of  1524  mm  long,  12.7  mm  diameter  incident  and  transmitted  bars  of  6061-T6  aluminum.  The 
striker  is  610  mm  long  and  made  of  the  same  material  as  the  other  bars.  The  samples,  whieh 
were  nominally  8  mm  diameter  by  3.5  mm  thiek  and  5  mm  diameter  by  2.5  mm  thick,  are 
positioned  between  the  incident  and  transmitted  bars.  The  bar  faees  were  lightly  lubricated  with 
grease  to  reduee  friction.  A  complete  description  of  this  testing  system  ean  be  found  in 
Reference  12. 

In  the  quasi-static  experiments,  the  elastic  modulus  was  determined  by  fitting  a  straight  line  to 
the  initial,  linear  part  of  the  stress-strain  eurve.  The  yield  stress,  in  both  quasi-static  and  dynamic 
experiments,  was  determined  by  fitting  a  second  order  polynomial  to  the  yield  region  of  the 
stress-strain  eurve  and  taking  the  derivative  to  determine  the  maximum. 

RESULTS  AND  DISCUSSION 

The  stess-strain  curves  for  each  material  are  shown  in  Figure  2  and  the  yield  stress  at  each  strain 
rate,  along  with  the  elastic  modulus  is  given  in  Table  2.  It  ean  be  seen  that  for  any  individual 
composite  the  stress  increases  with  strain  rate.  However,  it  is  more  difficult  to  determine  the 
effeet  of  the  volume  percent  of  aluminum  or  niekel  and  the  aluminum  partiele  size  purely  by 
observation. 


Table  2:  Elastic  modulus  and  yield  stress  for  particulate  composites,  MNML-1  through 
MNML-8,  with  varying  aluminum  and  nickel  volume  percents  and  aluminum  particle  size 


A1  Particle  Size 

Vf  A1 

VfNi 

E 

<7vs  (MPa)  1 

(urn) 

(GPa) 

£  =  0.0005  /s 

£  =  940  ±  150  /s 

£  =  4620  ±810/s 

MNML-l 

50 

40 

10 

7.0  ±  1.3 

91.0  ±0.6 

186  ±2 

210±2 

MNML-2 

5 

40 

10 

7.2  ±0.5 

97.7  ±3.7 

194  ±  1 

210±5 

MNML-3 

50 

20 

10 

5.2  ±0.6 

86.9  ±  2.2 

191  ±2 

213±3 

MNML-4 

5 

20 

10 

5.3  ±0.8 

92.1  ±0.4 

190  ±3 

208  ±2 

MNML-5 

50 

40 

0 

5.4  ±  1.4 

88.2  ±  1.9 

174  ±3 

188  ±5 

MNML-6 

5 

40 

0 

5.6  ±  1.7 

92.8  ±  1.0 

184  ±2 

203  ±3 

MNML-7 

50 
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Figure  2.  Stress-strain  charts  for  (a)  MNML-1,  (b)  MNML-2,  (c)  MNML-3,  (d)  MNML-4, 
(e)  MNML-5,  (f)  MNML-6,  (g)  MNML-7,  and  (h)  MNML-8 
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Blocks  of  material  mixed 
and  cast  in  random  order 
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and  tested  to  determine 
compressive  strength 


Figure  3.  Schematic  of  experimental  testing  and  analysis  as  a  split-plot  design 


The  materials  in  this  study  were  designed  using  a  factorial  design  of  experiments.  To 
accomplish  this,  two  levels  were  chosen  for  each  input  variable,  or  factor.  Each  possible 
combination  of  the  factors  was  then  prepared  as  an  MNML  particulate  composite.  The 
composites  were  prepared  in  random  order,  and  each  composite  was  cast  as  a  single  large  block 
of  material.  From  this  block,  several  samples  were  machined,  and,  subsequently,  approximately 
5  samples  of  each  material  was  tested  at  a  given  strain  rate,  as  shown  in  Figure  3.  This  is  not  a 
completely  randomized  factorial  design,  but  rather  a  split  plot  design  [14,15],  where  the  whole 
plot  factors  are  the  aluminum  particle  size  and  volume  percents  of  nickel  and  aluminum,  and  the 
subplot  factor  is  the  strain  rate.  This  design  allows  us  to  have  multiple  replicates  of  the  strain 
rate  factor.  However,  there  is  only  one  replicate  of  the  factors  in  the  whole  plot,  which  does  not 
allow  for  an  estimate  of  pure  error.  It  is  anticipated  that  similar  results  would  be  achieved  with  a 
completely  randomized  factorial  design,  versus  the  split  plot  design.  However,  fabrication  of 
two  smaller  blocks  of  material,  as  opposed  to  one  large  block,  would  have  permitted  replication 
of  the  whole  plot  and  provided  an  estimate  of  pure  error. 


The  elastic  modulus,  since  it  was  only  measured  at  one  strain  rate,  was  only  analyzed  as  a 
factorial  design.  The  ANOVA  table  of  results  is  given  as  Table  3,  where  a  significant  term  is 
defined  as  having  a  p-value  less  than  0.5.  The  adjusted  for  the  model  is  0.7291  indicating  that 
the  model  accountfs  for  nearly  73%  of  the  total  variability.  Interestingly,  only  the  volume 
fraction  of  nickel  and  aluminum  were  significant  factors.  The  aluminum  particle  size  and  any 
interactions  did  not  play  a  role  in  characterizing  the  elastic  modulus.  Figure  4  shows  the  effect 
of  increasing  the  aluminum  and  nickel  volume  percents,  which  is  to  increase  the  elastic  modulus. 
The  model  for  elastic  modulus,  in  terms  of  actual  factors,  is  given  by 

E  =  2.925  +  0.07*A  + 0.1 15*C  (2) 


where  A  is  the  volume  percent  of  aluminum,  C  is  the  volume  percent  of  nickel,  and  E  is  the 
elastic  modulus  in  GPa.  The  large  error  bars  on  the  effects  could  be  reduced  with  replicates  of 
the  materials  tested. 
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Table  3.  ANOVA  table  for  analysis  of  elastic  modulus 
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Figure  4. 


(a)  (b) 

Effect  of  (a)  aluminum  volume  percent  and  (b)  nickel  volume  percent  on  elastic 
modulus  in  aluminum-nickel  particulate  composites 


The  yield  stress  of  the  epoxy  binder  used  in  these  composites  has  been  found  to  demonstrate  a 
bi-linear  behavior  [16],  with  the  change  in  behavior  occuring  at  a  strain  rate  of  approximately 
100  /s.  Given  the  known  behavior  of  the  epoxy,  the  strain  rate  dependence  in  these  particulate 
composites  was  divided  into  low  strain  rate  and  high  strain  rate  regimes.  With  the  data  currently 
available  in  this  study,  the  effect  of  high  strain  rates,  i.e.  SHPB  rates,  can  be  investigated. 
However,  since  only  one  quasi-static  strain  rate  has  been  investigated  to  date,  the  effect  of  strain 
rate  at  low  rates  cannot  currently  be  analyzed.  Additional  testing  is  underway  to  determine  the 
effect  of  strain  rate  at  low  rates,  but  it  will  not  be  addressed  any  further  in  this  paper. 


The  analysis  for  the  yield  stress  at  high  strain  rates  is  slightly  more  complicated,  in  that  the  strain 
rate  is  analyzed  as  a  subplot  factor  and  the  remaining  factors  are  analyzed  as  whole  plot  factors. 
In  order  to  accomplish  this,  first  the  whole  plot  factors  are  analyzed,  neglecting  the  subplot  factor 
and  any  interactions  with  subplot  factor.  This  is  a  full  factorial  design  with  no  replicates.  Table 
4  gives  the  ANOVA  table  of  results  indicating  that  all  three  primary  factors  as  well  as  their 
interactions,  excluding  AC,  are  contributing  to  the  yield  stress  in  the  material. 


The  second  step  to  the  split  plot  analysis  is  to  analyze  the  split  plot  factors,  in  this  case,  strain 
rate.  In  this  analysis,  there  are  five  replicates,  since  each  material  was  measured  five  times.  This 
analysis  looks  at  the  primary  factor  of  strain  rate,  as  well  as  the  interactions  of  strain  rate  with  the 
factors  analyzed  as  part  of  the  whole  plot.  The  significant  factors  are  given  in  the  ANOVA  table 


for  this  analysis,  Table  5,  where  larger  F  vales,  and  resultingly  smaller  p  values,  indicate  higher 
levels  of  signficance  and  greater  contribution  to  the  model.  As  expected,  strain  rate  is  a 
significant  factor.  Interestingly,  strain  rate  interacts  with  the  volume  percent  of  aluminum  and 
nickel.  Also,  one  three  way  interaction  was  found  to  be  significant.  Since  the  strain  rate 
experiments  are  replicated,  the  residual  is  composed  of  both  lack  of  fit,  i.e.  the  error  due  to  the 
model  fitting,  and  pure  error,  which  is  the  error  associated  with  the  replicates.  It  can  be  seen 
that  the  lack  of  fit  error  is  not  significant,  indicating  that  the  model  sufficiently  describes  the 
available  data. 


Table  4.  ANOVA  table  for  analysis  of  yield  stress  whole  plot  factors,  where  A  is  the 

volume  percent  of  aluminum,  B  is  the  aluminum  particle  size,  C  is  the  volume 
percent  of  nickel,  and  terms  like  AB  indicate  interactions  between  to  of  the 
primary  factors 


Source 

Sum  of  Squares 

df 

Mean  Square 

F  value 

p-value 

Prob  >  F 

Model 

870.39 

5 

174.08 

64.31 

0.0154  significant 

A 

78.05 

1 

78.05 

28.83 

0.0330 

B 

75.78 

1 

75.78 

28.00 

0.0339 

C 

546.42 

1 

546.42 

201.88 

0.0049 

AB 

104.94 

1 

104.94 

38.77 

0.0248 

BC 

92.77 

1 

92.77 

34.27 

0.0280 

Residual 

5.41 

2 

2.71 

Cor  Total 

875.80 

7 

Now  that  the  analysis  on  the  whole  plot  and  split  plot  are  complete,  the  signifieant  factors  can  be 
combined  to  present  the  model  for  yield  stress  in  the  high  strain  rate  region, 

Yield  Stress  =  175.92039  +  0.26605  *  A_%A1  +  0.13312  *  B_AlSize  +  0.45416  *  C_%Ni 

+  5.88818E-003  *  D_strainrate  -  0.01 1 16  *  A_%A1  *  B_AlSize  -  6.12929E-005  (3) 

*  A_%A1  *  D_strainrate  +  7.87685E-003  *  B_AlSize  *  C_%Ni  +  4.70373E-005 

*  C_%Ni  *  D_strainrate  +  4.681 19E-006  *  B_AlSize  *  C_%Ni  *  D_strainrate, 

in  actual  units  for  the  factors  with  yield  stress  in  MPa.  Additional  tests  eould  be  conducted,  for 
example  fabricating  a  centerpoint  material,  in  order  to  validate  the  model. 

Figure  5  (a-d)  shows  the  effeets  of  the  primary  faetors  on  the  yield  stress,  with  the  remaining 
faetors  set  to  their  eenterpoint  values.  These  graphs  show  the  general  trends  for  the  factors,  but 
care  should  be  taken  sinee  all  the  faetors  were  involved  in  interaetions  and  the  graphs  are 
presented  at  levels  not  aetually  tested.  If  there  is  curvature  in  any  of  these  faetors,  it  will  affeet 
these  graphs.  The  aluminum  partiele  size,  the  volume  fraetion  of  niekel,  and  the  strain  rate  affect 
the  yield  stress  and  elastie  modulus  mainly  in  expeeted  ways,  based  on  the  load-sharing 
arguments  of  eomposite  theory.  As  the  niekel  volume  fraetion  increases  the  yield  strength 
inereases,  and  so  does  the  elastie  modulus.  The  yield  strength  also  increases  with  increasing 
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strain  rate.  Increasing  the  aluminum  volume  fraction  also  increases  the  elastic  modulus, 
according  to  Eqn.  (1),  albeit  with  a  smaller  factor  than  for  the  nickel.  This  is  consistent  with  the 
stiffer  elastic  modulus  of  the  nickel  compared  with  the  aluminum.  However,  as  the  aluminum 
volume  fraction  and/or  particle  size  is  increased,  the  yield  strength  decreases,  which  is  counter¬ 
intuitive.  One  possible  explanation  is  due  to  differences  in  the  strength  and/or  stiffness  of  the 
aluminum  particles  as  a  function  of  size,  due  to  the  presence  of  oxide  on  the  surface  of  the 
particles.  Another  explanation  is  that  this  could  be  due  to  chemical  interactions  between  the 
epoxy  and  the  aluminum,  however,  additional  experiments  (for  example  fabricating  a  centerpoint 
material)  would  be  needed  in  order  to  validate  this  and  any  alternative  hypotheses.  Figure  6  (a-d) 
shows  the  two  factor  interactions.  In  these  graphs,  one  factor  is  shown  on  the  x-axis  and  one 
factor  is  shown  as  two  different  colored  lines.  The  volume  fraction  of  aluminum  appears  to  have 
almost  no  effect,  within  the  error  bars,  at  small  particle  sizes  and  the  effect  of  decreasing  the 
yield  stress,  as  discussed  above,  at  the  large  particle  sizes.  Strain  rate  interacts  with  the  volume 
fraction  of  aluminum  and  nickel,  Figure  6  (b)  and  (d).  The  volume  fraction  of  nickel  interacts 
with  the  aluminum  particle  size.  This  may  indicate  that  the  addition  of  only  a  small  amount  of 
nickel  may  have  a  strong  effect  on  the  yield  strength. 

SUMMARY 

The  elastic  modulus  in  these  aluminum-nickel  particulate  composites,  was  found  to  depend  on 
the  volume  percent  of  both  the  nickel  and  the  aluminum  but  to  be  insenstive  to  aluminum  particle 
size,  within  the  bounds  of  the  factors  investigated.  The  yield  stress  at  high  strain  rates  was  found 
to  have  a  complex  dependance  on  all  of  the  factors  investigated  -  aluminum  and  nickel  volume 
percent,  aluminum  particle  size,  and  strain  rate.  If  a  desired  yield  strength  and  elastic  modulus 
are  known,  the  equations  developed  in  this  analysis,  i.e.  Equations  (1)  and  (2),  can  be  used  to 
simultaneously  optimize  the  properties  to  the  desired  amounts.  Additionally,  this  first  level 
analysis  could  be  used  to  further  refine  the  problem  with  the  addition  of  samples  to  determine  the 
possible  second  order  effects  of  the  factors. 

Further  testing  is  needed  at  quasi-static  strain  rates  in  order  to  determine  the  strain  rate 
dependance.  Additionally,  in  order  to  improve  the  factorial  design  of  experiments,  a  second 
block  of  materials  would  be  recommended.  In  this  block,  one  to  two  materials  would  be  repeated 
in  order  to  understand  the  block-to-block  variation  and  one  additional  material,  a  centerpoint, 
would  be  replicated  three  times  in  order  to  get  an  estimate  of  pure  error  as  well  as  an  estimate  of 
curvature.  This  second  block  would  require  the  manufacture  of  less  material  than  replicating  the 
entire  8  original  materials.  The  addition  of  a  centerpoint  would  also  help  to  elucidate  the  effect 
of  nickel  addition  on  the  composites,  as  the  current  design  has  either  no  nickel  or  some  nickel, 
which  is  different  than  the  change  in  volume  percent  of  aluminum  from  20%  to  40%. 
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Table  5.  ANOVA  table  for  analysis  of  yield  stress  split  plot  factor,  where  A,  B,  and  C  have 
the  same  meaning  as  Table  4,  and  D  is  the  strain  rate 


Source 

Sum  of  Squares 

df 

Mean  Square 

F  value 

p-value 

Prob  >  F 

Model 

6158.70 

4 

1539.68 

200.92 

<0.0001  significant 

D 

5984.46 

1 

5984.46 

780.96 

<0.0001 

AD 

72.81 

1 

72.81 

9.50 

0.0030 

CD 

132.77 

1 

132.77 

17.33 

<0.0001 

BCD 

127.93 

1 

127.93 

16.69 

0.001 

Residual 

490.43 

64 

7.66 

Lack  of  Fit 

62.48 

4 

15.62 

2.19 

0.0809  not  significant 

Pure  Error 

427.95 

60 

7.13 

Cor  Total 

6649.13 
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Figure  5.  Effect  of  (a)  volume  percent  of  aluminum,  (b)  aluminum  particle  size,  (c)  volume 
percent  of  nickel,  and  (d)  strain  rate  on  the  yield  stress  of  aluminum-nickel 
particulate  composites  under  high  strain  rate  loading 
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Figure  6.  Interactions  between  (a)  volume  percent  of  aluminum  and  aluminum  particle 
size,  (b)  volume  percent  of  aluminum  and  strain  rate,  (c)  volume  percent  of 
nickel  and  aluminum  particle  size,  and  (d)  volume  percent  of  nickel  and  strain 
rate  on  the  yield  stress  of  aluminum-nickel  particulate  composites  under  high 
strain  rate  loading 
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Abstract.  Multi-phase  particulate  composites  consist  of  individual  particles  of  more  than 
one  material  dispersed  throughout  and  held  together  by  a  polymer  binder.  The  mechanical 
and  physical  properties  of  the  composite  depend  on  the  properties  of  the  individual 
components;  their  loading  density;  the  shape  and  size  of  the  particles;  the  interfacial 
adhesion;  residual  stresses;  and  matrix  porosity.  These  multi-phase  particulate  composites 
systems,  particularly  those  with  high  fill  densities,  have  not  typically  been  studied  to 
determine  the  effects  of  microstructural  features  on  properties.  In  this  paper,  we  present  our 
investigation  of  the  influence  of  particle  size  and  dispersion  on  the  static  and  dynamic 
mechanical  response  of  these  multi-phase  (n>2)  polymer-metal  composites.  The  low  and 
high  strain  rate  compressive  strengths  are  determined  using  an  MTS  load  frame  and  a  split 
Hopkinson  pressure  bar,  respectively,  and  the  elastic  properties  were  studied  using 
dynamic  mechanical  analysis.  The  results  are  analyzed  using  a  factorial  design  of 
experiments  to  determine  the  effect  of  aluminum  and  nickel  volume  percent  and  aluminum 
particle  size  on  the  compressive  strength  as  a  function  of  strain  rate. 


1.  INTRODUCTION 

Particulate  composite  materials  composed  of  one  or  more  varieties  of  particles  in  a  polymer 
binder  are  widely  used  in  military  and  civilian  applications.  They  can  be  tailored  for  desired 
mechanical  properties  with  appropriate  choices  of  materials,  particle  sizes  and  loading  densities. 
Several  studies  on  similar  epoxy-based  composites  have  been  reported  and  have  shown  that 
particle  size  [1,  2],  shape  [3],  and  concentration  [4]and  properties  of  the  constituents  can  affect 
the  properties  of  particulate  composites.  In  composites  of  AI2O3  particles  in  epoxy  (Epon 
828/Z),  increasing  the  particle  concentration  and  decreasing  the  particle  size  is  found  to  increase 
the  strength  corresponding  to  4%  plastic  strain  [5].  A  study  of  aluminum  fdled  epoxy 
(DGEBA/MTHPA)  composites  has  found  that  a  small  amount  of  fdler  (~  5  vol.%)  increases  the 
compressive  yield  stress,  but  additional  amounts  of  fdler  decrease  the  compressive  yield  stress 
[6].  However,  tests  on  epoxy  (DOW  DER  331/bisphenol-A)  found  that  increasing  the  volume 
percent  of  glass  bead  fdler  increased  the  yield  stress  and  fracture  toughness  of  the  material  [7,  8]. 
In  a  study  on  a  similar  material,  decreasing  the  aluminum  particle  size  from  micro  to  nano 
resulted  in  increased  epoxy  crosslink  density  and  subsequently  increased  static  and  dynamic 
strength  [1].  In  this  study,  a  factorial  design  of  experiments  is  used  to  examine  the  effect  of 
aluminum  particle  size  and  aluminum  and  nickel  volume  percents  in  epoxy-based  Al-Ni 
particulate  composites. 
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2.  EXPERIMENTAL  SET-UP 


Composites  of  aluminum  and  nickel  powders  in  an  epoxy  binder  were  prepared  by  casting.  The 
aluminum  particle  size  was  varied  between  5  (Valimet,  H5  aluminum)  and  50  (Valimet,  H50 
aluminum)  pm.  The  H5  aluminum  was  found  to  have  an  average  particle  size  of  5.4  pm  with 
spherical  smooth  particle  morphology.  The  H50  aluminum  particles  were  also  smooth  and 
nominally  spherical  with  an  average  particle  size  of  51.9  pm.  The  nickel  particles  (Micron 
Metals)  had  rougher  surface  texture  and  more  irregular  shape,  with  an  average  particle  size  of 
47.5  pm.  However,  there  was  also  a  small  fraction  of  particles  in  the  nickel  powder  that  had  an 
average  particle  size  of  97.4  pm. 

Dynamic  mechanical  analysis  was  performed  using  a  TA  Q800  in  a  single  cantilever 
configuration  at  1,  10,  and  100  Hz  over  temperatures  from  148  K  to  473  K. 

Compression  experiments  at  quasi-static  strain  rates  (approximately  IxlO'"^  and  1x10'^  s'') 
were  conducted  with  an  MTS  810  testing  system  with  a  100  KN  test  frame  and  a  constant 
crosshead  displacement  rate.  Care  was  taken  to  center  the  samples  on  the  platens  prior  to  testing. 
A  thin  film  of  Boron  Nitride  (BN)  with  a  layer  of  Molybdenum  disilicide  (MOSi2)  on  top  was 
used  to  lubricate  the  surfaces  of  the  platen  in  contact  with  the  test  specimen.  In  addition  to  the 
MTS  system  recording  the  loads  and  displacement  of  the  frame  an  interfacing  software  (VIC 
Gauge  2.0  from  Correlated  Solutions  Inc.)  reads  input  voltages  for  both  the  load  and 
displacement.  This  software  interfaces  with  a  video  system,  which  allows  the  user  to  place 
virtual  displacement  gages  on  the  specimen  that  are  tracked  as  testing  takes  place.  A  high 
contrasting  boundary  or  point  is  required  for  tracking,  and  a  black  marker  was  used  to  draw 
fiducial  marks  on  the  specimens.  The  strain  is  measured  directly  from  the  specimen  itself  rather 
than  from  the  MTS  load  frame.  Multiple  virtual  displacement  gages  were  used  for  comparison 
and  to  enable  the  test  to  continue  in  the  event  that  one  gage  failed. 

Compression  experiments  at  intermediate  strain  rates  (approximately  1x10^  and  5x10^  s”') 
were  conducted  using  a  split  Hopkinson  pressure  bar  (SHPB)  [10]  system.  The  experiments 
were  conducted  using  the  SHPB  system  located  at  AFRL/RWME,  Eglin  AFB,  FL,  which  is 
comprised  of  1524  mm  long,  12.7  mm  diameter  incident  and  transmitted  bars  of  6061-T6 
aluminum.  The  striker  is  610  mm  long  and  made  of  the  same  material  as  the  other  bars.  The 
samples,  which  were  nominally  8  mm  diameter  by  3.5  mm  thick  and  5  mm  diameter  by  2.5  mm 
thick,  are  positioned  between  the  incident  and  transmitted  bars.  The  bar  faces  were  lightly 
lubricated  with  grease  to  reduce  friction.  A  complete  description  of  this  testing  system  can  be 
found  in  Reference  9. 

The  yield  stress,  in  both  quasi-static  and  dynamic  experiments,  was  determined  by  fitting  a 
second  order  polynomial  to  the  yield  region  of  the  stress-strain  curve  and  taking  the  derivative  to 
determine  the  maximum.  For  the  quasi-static  experiments,  the  yield  stress  reported  is  an  average 
of  the  values  from  all  of  the  virtual  strain  gauges. 


3.  RESULTS  AND  DISCUSSION 

The  composite  materials  used  in  this  study  were  prepared  using  a  factorial  design  of 
experiments,  with  the  variables  being  the  aluminum  particle  size  (5  or  50  pm),  the  volume 
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percent  of  aluminum  (20  or  40  vol.%),  and  the  addition  of  nickel  (10  vol.%).  The  full  factorial 
design  of  experiments  is  presented  in  Table  1.  Two  levels  were  chosen  for  each  input  variable, 
or  factor.  Each  possible  combination  of  the  factors  was  then  prepared  as  a  particulate  composite, 
labeled  MNML-1  through  MNML-8.  The  composites  were  prepared  in  random  order,  and  each 
composite  was  cast  as  a  single  large  block  of  material.  From  this  block,  several  samples  were 
machined,  and,  subsequently,  approximately  3-5  samples  of  each  material  were  tested  at  a  given 
strain  rate.  This  is  not  a  completely  randomized  factorial  design,  but  rather  a  split  plot  design 
[11,12],  where  the  whole  plot  factors  are  the  aluminum  particle  size  and  volume  percents  of 
nickel  and  aluminum,  and  the  subplot  factor  is  the  strain  rate.  This  design  allows  us  to  have 
multiple  replicates  of  the  strain  rate  factor.  In  contrast  to  the  subplot  factor,  there  is  only  one 
replicate  of  the  whole  plot  factors,  which  does  not  allow  for  an  estimate  of  pure  error.  It  is 
anticipated  that  similar  results  would  be  achieved  with  a  completely  randomized  factorial  design, 
versus  the  split  plot  design.  However,  fabrication  of  two  smaller  blocks  of  material,  as  opposed 
to  one  large  block,  would  have  permitted  replication  of  the  whole  plot  and  provided  an  estimate 
of  pure  error. 


Table  1.  Yield  stress  for  particulate  composites,  MNML-1  through  MNML-8,  with  varying 
aluminum  and  nickel  volume  percent  (V%)  and  aluminum  particle  size,  where 
_ the  headers  for  the  yield  stress  column  indicate  the  strain  rate _ 


A1  Particle 

Ovs  (MPa) 

Size 

(pm) 

V%  A1 

V  /o 

Ni 

9a:10-^s-i 

950S-1 

4600s- 

MNML-1 

50 

40 

0.1 

89.6  ±  0.4 

102  ±  1 

186  ±2 

210±2 

MNML-2 

5 

40 

0.1 

94.6  ±  2.7 

104  ±2 

194  ±  1 

210±5 

MNML-3 

50 

20 

0.1 

86.5  ±3.5 

90.8  ±  1.1 

191  ±2 

213±3 

MNML-4 

5 

20 

0.1 

86.0  ±0.6 

98.6  ±2.5 

190  ±3 

208  ±2 

MNML-5 

50 

40 

0 

82.6  ±2.2 

94.4  ±2.7 

174  ±3 

188  ±5 

MNML-6 

5 

40 

0 

92.0  ±0.8 

100  ±3 

184  ±2 

203  ±3 

MNML-7 

50 

20 

0 

83.4  ±  1.4 

97.5  ±  0.6 

183  ±2 

198  ±2 

MNML-8 

5 

20 

0 

85.6  ±  1.1 

99.2  ±3.5 

187  ±2 

202  ±5 
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Figure  1.  (a)  Stress-strain  curve  for  MNML-5  and  (b)  peak  stress  as  a  function  of  strain 
rate  for  MNML-5 


An  example  stress-strain  eurve  for  MNML-5  is  shown  in  Figure  1  (a),  and  the  yield  stress  at 
each  strain  rate  is  given  in  Table  1  with  a  representative  curve  given  in  Figure  1  (b).  It  can  be 
seen  that  for  any  individual  composite  the  peak  stress  increases  with  strain  rate.  However,  it  is 
more  difficult  to  determine  the  effect  of  the  volume  percent  of  aluminum  or  nickel  and  the 
aluminum  particle  size  purely  by  observation. 

The  yield  stress  of  the  epoxy  binder  used  in  these  composites  has  been  found  to  demonstrate 
a  bi-linear  behavior  [13],  with  the  change  in  behavior  occuring  at  a  strain  rate  of  approximately 
200  s'^  as  determined  by  yield  stress  measurements  in  compression.  Given  the  known  behavior 
of  the  epoxy,  the  strain  rate  dependence  in  these  particulate  composites  was  divided  into  low 
strain  rate  and  high  strain  rate  regimes.  In  order  to  verify  that  this  separation  is  valid,  dynamic 
mechanical  analysis  was  conducted  on  all  samples.  The  elastic  modulus  at  room  temperature  for 
strain  rates  from  lO'"^  -  10^  was  determined  using  the  decompose/shift/reconstruct  (DSR)  method 
described  by  Mulliken  and  Boyce  [14].  A  straight  line  was  then  fitted  to  the  low  and  high  strain 
rate  portions  of  the  curve  and  the  intersection  of  these  two  lines  was  taken  as  the  strain  rate  at 
which  the  beta  transition,  which  in  many  glassy  polymers  is  a  low  temperature  transition 
corresponding  to  side  chain  motion,  moves  to  room  temperature.  The  transition  strain  rate  was 
analyzed  as  a  reponse  to  the  whole  plot  factors,  and  the  model  was  determined  to  be 
insignificant,  i.e.  the  values  are  the  same  within  the  statistical  error.  The  average  transition  strain 
rate  for  MNML-1  through  -8  is  10.5  s'*,  where  the  same  transition  determined  for  epoxy  is  10  s'*. 
Although  the  lack  of  significance  could  be  due  to  large  system  noise  or  uncontrolled  variables 
that  are  having  effect,  since  the  DMA  test  is  looking  primarily  at  the  properties  of  the  polymer 
binder,  this  analysis  suggests  that  the  epoxy  transitions  are  not  being  changed  by  the  factors. 
This  analysis  verifies  the  validity  of  separating  the  analysis  into  a  low  and  high  strain  rate  regime 
and  defines  the  appropriate  areas  for  both  regimes.  Figure  2  (b)  shows  the  transition  strain  rate 
along  with  the  peak  stress  for  MNML-5. 

The  analysis  for  the  yield  stress  in  both  strain  rate  regimes  is  slightly  complicated,  in  that 
the  strain  rate  is  analyzed  as  a  subplot  factor  and  the  remaining  factors  are  analyzed  as  whole 
plot  factors.  In  order  to  accomplish  this,  first  the  whole  plot  factors  are  analyzed  for  each  strain 
rate  regime,  neglecting  the  subplot  factor  and  any  interactions  with  subplot  factor.  This  is  a  full 
factorial  design  with  no  replicates.  Tables  2  and  3  gives  the  ANOVA  (analysis  of  variance)  table 
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of  results  for  the  low  strain  rate  regime  and  the  high  strain  rate  regime,  respectively.  For  the  low 
strain  rate  regime,  4  factors  and  one  interaction  are  contributing  to  the  yield  stress  of  the 
material,  indicated  by  a  p-value  less  than  0.05.  Although  volume  percent  of  nickel  has  a  p-value 
greater  than  0.05  indicating  that  it  is  not  significant  to  the  model,  it  is  included  for  hierarchy, 
since  the  interaction  AC  is  significant.  It  is  unlikely  that  an  interaction  would  be  significant  if 
the  main  effect  is  not  playing  a  role.  For  the  high  strain  rate  regime,  all  three  primary  factors  as 
well  as  their  interactions,  excluding  AC,  are  contributing  to  the  yield  stress  in  the  material. 

The  second  step  to  the  split  plot  analysis  is  to  analyze  the  split  plot  factors,  in  this  case, 
strain  rate.  In  this  analysis,  there  are  3  -  5  replicates,  since  each  material  was  measured  multiple 
times.  This  analysis  looks  at  the  primary  factor  of  strain  rate,  as  well  as  the  interactions  of  strain 
rate  with  the  factors  analyzed  as  part  of  the  whole  plot.  The  significant  factors  are  given  in  the 
ANOVA  tables  for  the  low  and  high  strain  rate.  Tables  4  and  5,  respectively.  As  expected,  strain 
rate  is  a  significant  factor.  Interestingly,  at  low  strain  rates,  the  strain  rate  does  not  appear  to 
interact  with  any  of  the  other  factors,  i.e.  the  interactions  all  had  p-values  over  0.05  making  them 
insignificant  for  the  model.  However,  at  high  strain  rates,  the  strain  rate  interacts  with  the 
volume  percent  of  aluminum  and  nickel.  Also,  one  three  way  interaction  was  found  to  be 
significant.  Since  the  strain  rate  experiments  are  replicated,  the  residual  is  composed  of  both 
lack  of  fit,  i.e.  the  error  due  to  the  model  fitting,  and  pure  error,  which  is  the  error  associated  with 
the  replicates.  It  can  be  seen  that  the  lack  of  fit  error  is  not  significant,  indicating  that  the  model 
sufficiently  describes  the  available  data. 

Table  2.  ANOVA  table  for  analysis  of  yield  stress  whole  plot  factors  at  low  strain  rates, 

where  A  is  the  volume  percent  of  aluminum,  B  is  the  aluminum  particle  size,  C  is 
the  volume  percent  of  nickel,  and  terms  like  AB  indicate  interactions  between  the 
primary  factors _ ^ _ _ _ 


Source 

Sum  of  Squares 

df 

Mean  Square 

F  value 

p-value 

Prob  >  F 

Model 

141.32 

4 

35.33 

16.60 

0.0219  significant 

A 

46.50 

1 

46.50 

21.84 

0.0185 

B 

60.86 

1 

60.86 

28.59 

0.0128 

C 

3.30 

1 

3.30 

1.55 

0.3017 

AC 

29.91 

1 

29.91 

14.05 

0.0332 

Residual 

6.39 

3 

2.13 

Total 

147.70 

7 
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Table  3.  ANOVA  table  for  analysis  of  yield  stress  whole  plot  factors  at  high  strain  rates, 
where  A  is  the  volume  percent  of  aluminum,  B  is  the  aluminum  particle  size,  C  is 
the  volume  percent  of  nickel,  and  terms  like  AB  indicate  interactions  between  the 
primary  factors _ ^ _ _ _ 


Source 

Sum  of  Squares 

df 

Mean  Square 

F  value 

p-value 

Prob  >  F 

Model 

870.39 

5 

174.08 

64.31 

0.0154  significant 

A 

78.05 

1 

78.05 

28.83 

0.0330 

B 

75.78 

1 

75.78 

28.00 

0.0339 

C 

546.42 

1 

546.42 

201.88 

0.0049 

AB 

104.94 

1 

104.94 

38.77 

0.0248 

BC 

92.77 

1 

92.77 

34.27 

0.0280 

Residual 

5.41 

2 

2.71 

Total 

875.80 

7 

Table  4.  ANOVA  table  for  analysis  of  yield  stress  split  plot  factor  at  low  strain  rates, 
where  A,  B,  and  C  have  the  same  meaning  as  Table  2,  and  D  is  the  strain  rate 


Source 


Sum  of  Squares 


df 


Mean  Square 


F  value 


p-value 
Prob  >  F 


Model 


1857.29 


464.32 


86.90 


< 


0.0001  significant 


D 


177 S. 98 


1773.98 


332.02 


<  0.0001 


CD 


56.61 


56.61 


10.60 


0.0020 


ABD 


48.24 


48.24 


9.03 


0.0041 


ACD 


Residual 


Lack  of  Fit 


Pure  Error 


Total 


28.71 


0.0246 


267.15 


34.14 


233.01 


2124.44 


Table  5.  ANOVA  table  for  analysis  of  yield  stress  split  plot  factor  at  high  strain  rates, 
where  A,  B,  and  C  have  the  same  meaning  as  Table  3,  and  D  is  the  strain  rate 


Source 


Sum  of  Squares 


df 


Mean  Square 


F  value 


p-value 
Prob  >  F 


Model 


6158.70 


1539.68 


200.92 


<0.0001  significant 


D 


5984.46 


5984.46 


780.96 


<0.0001 


AD 


72.81 


72.81 


9.50 


0.0030 


CD 


132.77 


132.77 


17.33 


<0.0001 


BCD 


Residual 


Lack  of  Fit 


Pure  Error 


Total 


127.93 


127.93 


16.69 


0.001 


490.43 


62.48 


427.95 


6649.13 
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Now  that  the  analysis  on  the  whole  plot  and  split  plot  are  complete,  the  significant  factors 
can  be  combined  to  present  the  experimental  design  model  for  yield  stress  in  the  low  strain  rate 
region, 

Yield  stress  =  84.67029  +  0.054922  *  A_%A1  -  0.079909  *  B_AlSize  -  0.34211  *  C_%Ni 
(1) 

+16828.15774  *  D_strainrate  +  0.020227  *  A_%A1  *  C_%Ni 

-  689.06577  *  C_%Ni  *  D_strainrate  -  0.38121  *  A_%A1  *  B_AlSize  * 

Dstrainrate 

+  11.40636  *  A_%A1  *  C_%Ni  *  D_strainrate, 
and  the  high  strain  rate  region. 

Yield  stress  =  175.92039  +  0.26605  *  A_%A1  +  0.13312  *  B_AlSize  +  0.45416  *  C_%Ni 
(2) 

+  5.888 18E-003  *  D_strainrate  -  0.01116  *  A_%A1  *  B_AlSize 

-  6.12929E-005  *  A_%A1  *  D_strainrate  +  7.87685E-003  *  B_AlSize  *  C_%Ni 

+  4.70373E-005  *  C_%Ni  *  D_strainrate 

+  4.681 19E-006  *  B  AlSize  *  C_%Ni  *  D  strainrate, 

in  actual  units  of  the  factors  with  yield  stress  in  MPa.  It  is  interesting  that,  although  all  of  the 
main  effects  were  found  to  be  significant  at  both  strain  rates,  the  interactions  were  not  found  to 
be  the  same.  Neglecting  the  interactions,  for  both  strain  rates,  the  aluminum  particle  size,  the 
volume  fraction  of  nickel,  and  the  strain  rate  affect  the  yield  stress  mainly  in  expected  ways, 
based  on  a  load  sharing  mechanism  between  the  much  stiffer  metal  particles  and  the  epoxy 
matrix.  At  low  strain  rates,  as  the  aluminum  and/or  nickel  volume  fraction  increases,  the  yield 
strength  increases.  Increasing  the  strain  rate  also  increases  the  yield  strength  of  the  material  -  a 
similar  effect  to  that  which  is  observed  even  in  unreinforced  epoxies.  Increasing  the  aluminum 
particle  size  leads  to  a  reduction  in  yield  strength  -  this  could  be  due  to  differences  in  the 
strength  and/or  stiffness  of  the  aluminum  particles  as  a  function  of  size,  due  to  the  presence  of 
oxide  on  the  surface  of  the  particles.  Less  satisfactory,  however,  is  the  decrease  in  yield  strength 
with  increasing  aluminum  volume  fraction  at  high  strain  rates.  It  is  possible  that  this  could  be 
due  to  chemical  interactions  between  the  epoxy  and  the  aluminum,  however,  additional 
experiments  (for  example  fabricating  a  centerpoint  material)  would  be  needed  in  order  to 
validate  this  hypothesis. 

4.  SUMMARY 

The  yield  stress  at  low  and  high  strain  rates  was  found  to  have  a  complex  dependance  on  all 
of  the  factors  investigated  -  aluminum  and  nickel  volume  percent,  aluminum  particle  size,  and 
strain  rate.  If  a  desired  yield  strength  as  a  function  of  strain  rate  is  known,  the  equations 
developed  in  this  analysis,  i.e.  Equations  (I)  and  (2),  can  be  used  to  simultaneously  optimize  the 
properties  to  the  desired  amounts.  Additionally,  this  first  level  analysis  could  be  used  to  further 
refine  the  problem  with  the  addition  of  samples  to  determine  whether  the  factors  exhibit  second 
order  effects. 
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In  order  to  improve  the  factorial  design  of  experiments,  a  second  block  of  materials  would 
be  recommended.  In  this  block,  one  to  two  materials  would  be  repeated  in  order  to  understand 
the  block-to-block  variation  and  one  additional  material,  a  centerpoint,  would  be  replicated  three 
times  in  order  to  get  an  estimate  of  pure  error  as  well  as  an  estimate  of  curvature  in  the  whole 
plot  factors.  This  second  block  would  require  the  manufacture  of  less  material  than  replicating 
the  entire  8  original  materials.  The  addition  of  a  centerpoint  would  also  help  to  elucidate  the 
effect  of  nickel  addition  on  the  composites,  as  the  current  design  has  either  no  nickel  or  some 
nickel,  which  is  different  than  the  change  in  volume  percent  of  aluminum  from  20%  to  40%. 
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